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Abstract 

We present a Borel resummation method for the QCD corrections in inclusive, charmless, 
semileptonic B meson decay. The renormalon contributions are resummed to all orders by em- 
ploying a bilocal expansion of the Borel transform that accurately accounts for the first infrared 
renormalons in the Borel plane. The renormalons in the pole mass and the QCD expansion are 
resummed separately, and a precise relation is obtained between a properly defined pole mass and 
the MS mass. The inclusive decay rate is calculated to three loop order using an estimate of the 
yet unknown NNLO coefficient. 
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The Cabibbo-Kobayashi-Maskawa matrix element |Kib|, which is an important parameter 
in the study of CP violation in weak interactions, can be determined through the inclusive, 
charmless, semileptonic B decay. Presently, the largest theoretical uncertainty in the decay 
rate comes from the b quark mass and the QCD correction. It is well known that these two 
are not independent but intimately related. In the large b quark mass limit, the B meson 
decay rate is given by that of a free b quark with a correction quadratically suppressed in 
the b mass The decay rate in the on-shell scheme is given in the form 

TiB X„,,J = ^f^rnlf{a,)[l + 0(l/m^)] (1) 

with mb and /(os) denoting the b quark pole mass and the QCD correction to the free b 
quark decay, respectively. In perturbation theory /(os) can be expanded as 

oo 

/ = i+5:/„(eK(/ir+\ (2) 

n=0 

where ^^(/i) is the strong coupling constant and ^ = fi/m^, with mj^[= "^ms("^ms)] 
noting the MS mass. 

Because of the infrared (IR) renormalon the perturbative coefficients fn are expected 
to grow rapidly, and this would give rise to an O^Aqc-D/fnb) intrinsic uncertainty to the 
asymptotic series. On the other hand, this renormalon uncertainty in /(a^) is known to 
arise from the use of the pole mass in the decay rate. When the pole mass is expanded 
in the strong coupling it also has a renormalon caused uncertainty of O(Aqcd), and this 
uncertainty cancels that of /(os) 0, 0|. Once the pole mass is replaced by a renormalon 
free mass like the MS mass or one of the threshold masses 0, 0, 0, B the resulting /(a,) 
is in principle expected to give a well-behaved series. In practice, however, the use of the 
renormalon free masses are not completely satisfactory. For instance, with the MS mass the 
renormalon cancellation is not so obvious at low orders and the convergence for /(a^) is 
uite slow, possibly due to the large power of nib in the decay rate, and with the IS* mass 
one must assume that the 15* state of T is a perturbative system. 
In this paper we describe a more direct approach to the renormalon problem, in which 
the leading renormalons in the pole mass and /(a^) are Borel resummed separately. An 
essential ingredient of this program is the accurate calculation of the renormalon residue for 
the pole mass. The residue, which determines the normalization constant of the large order 
behavior, can be computed in perturbation theory and it turns out that the residue 

for the pole mass can be determined to good accuracy within a few percent error 11, ^ 13l |. 



The renormalon cancellation then allows the residue for /(os) to be determined to the same 
accuracy from the pole mass residue. The Borel transforms of the pole mass and / can be 
systematically expanded around the origin and the first renormalon location in the Borel 
plane. With the computed residues these two expansions can be interpolated to a bilocal 
expansion to obtain an accurate description of the Borel transforms in the region in-between 
the two expansion points. Since at the typical strong coupling in B decay the bulk of the 
Borel integral comes from this region of integration, the bilocal expansion can much improve 
the accuracy of the Bore resummation. By taking into account the renormalon singularity 
correctly this program essentially resums the renormalon caused large order behavior to all 
orders. 

With the renormalon singularity on the integration contour one has to worry about the 
ambiguities of O(Aj^) for the pole mass and 0{Aj^/mb) for / arising from the ambiguities 
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in choosing the integration contour that may either be on the upper or the lower half 
plane beginning at the renormalon singularities. It was shown in Refs. that these 

ambiguities cancel as long as one takes the same contour for the pole mass and 
Furthermore the authors argued that in general the pole mass, when consistently used, 
does not cause any ambiguity to physical observables in the heavy quark effective theory. 
Until recently, however, there was no procedural program implementing this observation. 
The main reason for this, perhaps, was that one did not know how to describe the Borel 
transforms beyond the immediate neighborhood of the origin in the Borel plane. Within 
the large /3o approximation 0] one can obtain the Borel transforms in the whole Borel 
plane, but the accuracy of the approximation is limited. What is necessary for practical 
applications is an accurate description of the Borel transforms in the region in-between the 
origin and the first renormalon singularity, since the contribution from the region far from 
the origin is exponentially suppressed. The bilocal expansion fills this gap. 

The program described above was successfully applied to the heavy quark static potential 
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which is known to suffer from a severe convergence problem with the ordinary perturba- 
tive expansion. The Borel resummed potential obtained from employing a bilocal expansion 
turns out to have an excellent convergence and is in remarkable agreement with the lattice 
calculations. 

Throughout the paper, unless stated otherwise, the perturbative expansions are assumed 
to be in the MS scheme with four active quark flavors {Nf = 4), and only the leading IR 
renormalon closest to the origin in the Borel plane is considered. 

The pole mass in terms of the MS mass mjj^ has the perturbative expansion 



rrib = 



1 + 



oo 

E 

n=0 



Pn(0«s(/i) 



n+1 



(3) 



The first three coefficients Pn are known at fi = rriy^ fl3, E, 1^ and their /i dependence 
can be obtained by the renormalization group invariance of the pole mass. The Borel 
resummation of this expansion can be written formally as 



mb[as{fx)] = 



ooitie 



-b/l3oasifJ.)~ 



Th{b,^) db 



where the Borel transform m{b,S,) has the perturbative expansion 

KM ' 



(^0 = E 



n=0 



(4) 



(5) 



with /3o denoting the one loop coefficient of the QCD {3 function [(5 = das/ dYn.^jji'^) = — a^(/3o + 
I3ias+(32a1 + ■ ■ ■)]■ The integration contour, which can be derived from the general argument 
of Borel resummation of a series with a same sign large order behavior can be either 
on the upper or the lower half plane along the positive real axis. The pol e mass is known to 
have an ambiguity from renormalon singularity proportional to 21| that has the weak 
coupling expansion 



^MS ^ f^Ctsifi)' 



1 + ria^(/i) + r2as{iif + 



(6) 



where 



2 ' 



ri 
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For the Borel integral to generate this ambiguity one can easily see that the Borel 
transform m{b,C,) must have the singularity 



C f 

Mb, = [l + ci(l - 2b) + C2(l - 2bf + ..]+ analytic part , 



where 

The 'analytic part' denotes any terms analytic around 6=1/2 with radius of convergence 
bounded by the next renormalon at 6 = 3/2 and Cm is a real constant. 

With this singularity in fh{b,^) the resummed pole mass can be written formally as 

rrib = mBR[as{fi)] ± iTm[as{fi)] , (10) 

where the 'Borel resummed (BR)' mass denotes the real part of the Borel integral in 
Eq. Q and denotes the imaginary part obtained with the integration contour on the 
upper half plane. Note that by definition the BR mass has the same perturbative expansion 
as the pole mass and is not a short distance mass. 

The residue Cm, which is not known exactly, can be calculated systematically in pertur- 
bation theory [ol. [lol|. The details of the calculation adopted here can be found in 11 1. To 



compute the residue we first expand using the known up to NNLO and the unknown 
the function R{b,^) = (1 — 26)''^^'^^m(6, ^) about the origin in the w plane defined by the 
following conformal mapping 



vi+^ - yi - 2&/3 

w = , (11) 

Vm + Jl - 26/3 



which gives to 0{w^) 

R[b{w), 1] = 0.42441 + 0.61198W + 0.25351^^ + [-108.99217 + 7.90052p3(l)]w;^ . (12) 
This allows us to estimate p^ i^l . From the pattern of the known lower order terms we may 



safely assume |108.99217 - 7.90052 p3(l)| < 1 , which leads to 

P3(l) = 13.796 ±0.127. (13) 

Now evaluating Eq. ()12|) at the renormalon location w = l/5 we obtain the residue 

Cm = R[b{w = 1/5),1] ^ 0.42441 + 0.12240 + 0.01014 ±0.00800 = 0.55695 + 0.00800, (14) 

where the error comes from the uncertainty in the estimate (|13|). 

The BR mass can be obtained accurately by an interpolation of the two expansions 
(0) and (jHI) of the Borel transform rfi(6, ^) (bilocal expansion [ll|). Since at the typical values 
of the strong coupling considered here the bulk of the contribution to the Borel integral comes 
from the region between the origin and the first IR renormalon, it is important to have the 
Borel transform in this region as accurately as possible from the known first perturbative 
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terms. As demonstrated in the heavy quark potential the accuracy of the Borel transform 
can be improved greatly through the bilocal expansion 

The bilocal expansion of the Borel transform for the pole mass is given by 



m{b,^) = lim mN,M(&,0 



N,ll^^\^^ n\ [PoJ (1-26)1+ 



l + ^Q(l-26r 



i=l 



(15) 



where g„ are to be determined by demanding Eq. (fTH|) reproduce Eq. ^ when expanded 
about 6 = 0. With the known first two coefficients Ci^2, and taking M = 2, we have the first 
three g„ as 

?o(0 = PO- Cm^{l + Ci + C2) , 

qiiO = PliO - '2CUPo[l - C2 + z/(l + ci + C2)] , 

92(0 = P2{0 - 4«/5o'[2 + ^^(3 + ci - C2) + z/2(l + ci + C2)] . (16) 

Since the two expansions about the origin and about the renormalon singularity are expected 
to be convergent on the disks, |6| < 1/2 and \b — l/2\ < 1, respectively, the interpolation 
(I15|) is expected to give a good description of the Borel transform in the region that includes 
the two expansion points. 

Taking the canonical value for the strong coupling as{mj^) = 0.22 and ^ = 1, and 
substituting the known rriN 2(&, 1)) (^ = 0; I5 2), into the Borel integral (H)), and performing 
the integration in the w-plane for the convenience of numerical integration, we obtain 

= m^{l + 0.15769 + 0.00409 - 0.00028 T 0.00014) , (17) 

where the error was obtained by varying Cm within the uncertainty in the computed residue 
in Eq. p4|l . Notice the remarkable convergence of the BR mass compared to the power 
series ® which at ^ = 1 gives rUb = m^{l + 0.09337 + 0.04550 + 0.03235) . This shows that 
the bilocal expansion improves the convergence by two orders of magnitude. 

We now turn to the resummation of /(a^). For this we first note that because of the large 
power {ml) in the pole mass in the decay rate the expansion Q may not be an optimal 
expansion for /(a^)- We can reorganize the expansion (j2)) by writing 

f[a,{fx)] = Fil + H[a,{fx)]) (18) 

where F{x) is an arbitrary function analytic at x = 1. Since H itself can be regarded as a 
coupling constant, with its relation to the MS coupling constant defined implicitly by Eq. 
fll8|) . we see that a particular choice of F is nothing but a renormalization scheme choice. 
Given an F the perturbative expansion for H 

00 

H=Y1 hniOc^sifiT"-' (19) 

71=0 

can be easily obtained using the definition ()18p and the series (j2I) for /(os). Instead of 
resumming the series directly we shall resum the expansion ()19|1 for H. 

The Borel resummation for H can be done in a similar fashion as in the pole mass resum- 
mation, assuming that the first IR renormalon for H is known. Then with the resummed 
pole mass and H we have 

mlf = {ruBR ± ^r^)'F(l + Hbk ± iTh) , (20) 
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where Hbr denotes the real part of the Borel integral 



H[as{fi)] 



1 



ooztie 



H{b,0 db 



(21) 



and Th denotes the imaginary part obtained with the contour on the upper half plane. The 
Borel transform H[b,^) has the perturbative expansion 



n=0 



'b^ 

Jo, 



(22) 



The structure of the IR renormalon for H can be easily obtained by the cancellation of 
the renormalons in ml and /(os). Expanding Eq. (^01) in powers of Tm,H we obtain 



;/ = ml^F{l + i^BR) ± t [5F(1 + HBK)Tm + F'{1 + ifBR)mBRr^,] + 0(r^.^) (23) 



The cancellation of the renormalons means vanishing of the ambiguous imaginary part, 
which gives 

5F[l + i7BR(as)] ^ , ^24) 



F'[l + i^BR(as)]"^BR(tts 



[a. 



The renormalon singularity for H{b,C,) can now be easily obtained by expanding the right- 
hand- side of this equation in small a^. 

At this point we choose a specific form for F{x). Considering the large power of the pole 
mass in the decay rate a plausible choice for F is 



F{x) 



l + -(x- 1) 
Q 



(25) 



with q a number not necessarily an integer. F{1 + H) is then essentially a power function 
with a rescaling of the normalization for H that is for the purpose of convenience only and 
does not have any physical effect. 

Using the perturbative expansions in Eqs. © and (|T!Hl of ttibr and i^BR (which have the 
identical expansions with mb and if, respectively) and obtained from the renormalon 
singularity (jHI), we have 



5Cm^r(— z/) sin(i/7r) 



1 + rias + r2{C)al + 



(26) 



where 



ri = ri -po + ho/q 
hiO = r2-piiO+PoiPo-ri) + [hi{^) + hori-poho]/q, 



(27) 



with ri^2 given in Eq. (jZj). For the Borel integral (PT|) to have this imaginary part the Borel 
transform H{b,^) must have the singularity 



H{b,0 



;i - 26) 



1 + ci(l -26) + 52(0(1 -26)' 



+ analytic part 



(28) 
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FIG. 1: The ratio of the computed renormalon residue to the expected one from the renormalon 
cancellation at NLO (solid) and at estimated NNLO (dashed) vs renormalization scale ^ = fi/rriy^. 

with 

C„K) = -5a.«. c, = J^, J,K) = j-^ML_. (29) 

Now the decay rate in the BR scheme is given by 

r(6 ^ X„,,J = 9M^rniM<^s)[l + 0{l/ml^)] , (30) 

where /br = F{1 + Hbr). Note that this is nothing but the decay rate in the on-shell 
scheme with Borel resummation, but now there is no theoretical ambiguity inherent in the 
formula in Eq. (P). Unlike /br, for instance, /(««) in Eq. ((T)) is not well defined due to the 
divergence of the expansion (j21). 

We now take the q ^ oo limit in Eqs. and ^T7\ . which is equivalent to putting 

/ = exp{H). What is nice with this limit is that the renormalon cancellation now occurs 
within the sum of the series for H and that of 51n(mfe). We think this is an attractive 
feature, considering the good cancellation of renormalons within the sum of the pole mass 
and the static potential in heavy quarkonium system 0, Also, we shall see that the q 
dependence of /(a^) becomes minimal in this limit. 

Before starting computation of -f^BR we first check how well the renormalon residue Ch{C,) 
computed using the known NLO calculation for H compares to the one in Eq. ()29|) expected 
from the renormalon cancellation. This exercise will allow us to estimate the unknown NNLO 
coefficient for /. The calculation for the residue Ch{0 is completely similar to that of the 
pole mass residue Cm- The ratio of the calculated NLO residue to the expected residue 
is plotted in Fig. [TJ The best matching between the calculated and the expected occurs 
at ^ ~ 0.34 with about 80% agreement. This scale is smaller than b quark mass, and 
this is in line with the expectation that the optimal expansion for /(og) would occur at 
a scale somewhat below the b mass. This also shows that the NNLO contribution to the 
residue should be significant. Now the unknown NNLO coefficient /2 can be estimated in 
the following way. We demand the NNLO residue at ^ = 0.34 fall within a conservative 20% 
of the expected value (|ratio — 1| < 0.2), which gives 

0.8 < -1.552-0.277/2(1) < 1.2. (31) 
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FIG. 2: /br at LO (dotted), NLO (dashed), and estimated NNLO (solid) at as(m,j^) = 0.22 vs 
renormalization scale ^ = /i/mj^. The shaded band shows the uncertainty of /br at NNLO. 

This provides us an estimate of the NNLO coefficient: 

/2(1) = -9.22 ±0.72. (32) 

This corresponds to / = [1 - 2.41a, - 21.300^ - (286 ± 22)0^] and mf / = m^[l + 4.25a, + 
26.78a^ + (160 ± 22)a^] where a, = as{'mj^)/7T. The latter shows that our estimate of the 
NNLO coefficient sits in the middle of the estimate 5^ = 125 based on a naive power growing 
behavior I23J and the Pade estimate 188 0|. 

The computation of -f^BR is completely parallel to that of ^br. We first combine the 
two expansions in Eqs. ()22|1 and (j^H|l for H into the bilocal expansion HN,M{b,C,) defined 
similarly to that of the pole mass in Eq. (jl5|) . The first two iJjv,2(&;0 = 0)1) can 
be obtained from the known NLO expansions (25| of /(a,), and the NNLO H2,2{b,C) from 
the estimated coefficient in Eq. (jH^ . Substituting these into the Borel integral (PT|) and 
performing the integration the Borel resummed -f^BR up to NNLO is obtained. The result for 
/br = exp (-^br) is in Fig. |21 The plots show that the scale dependence is significant at LO 
but becomes milder at NLO and small at NNLO for ^ > 0.5. Note that the NLO result is 
completely within the NNLO band arising from the error in the /2 estimate ()32|) . From the 
perspective of scale dependence the estimate in Eq. (jH^ with the lower bound in the error 
estimate is preferred because it defines the lower bound of the NNLO band which has a less 
scale dependence than the upper one. In the presence of this unphysical scale dependence 
we can take the principle of minimal sensitivity (PMS) j2^ values as the predictions of the 
NNLO resummations. This gives 

/bj^ = e^BR = 0.728 ± 0.011 , (33) 

where the error was taken from the difference between the PMS values of the central line and 
the lower bound of the NNLO band. Since an error estimate is always subjective, a more 
conservative estimate can be obtained, for example, by taking the variation of the NNLO 
estimates at ^ = 0.4 as the uncertainty, which gives /br ~ 0.728 ± 0.017. 

We now come back to the case of a finite value of q. In principle /br(q;s) should not 
depend on g, but with a finite order perturbation there is always scheme dependence, and a 
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FIG. 3: Renormalization scheme dependence of /br at NNLO (estimated). 



scheme with a less scheme dependence is obviously preferable. To see the q dependence we 
repeated the NNLO computation of /br at various values of g, and the result is in Fig. 01 
The scheme dependence at small q is quite significant, for example, the difference between 
g = 1 and g = 5 is larger than the NNLO contribution at ,^>0.25. For large q the scheme 
dependence is very mild and becomes minimal at g = oo. This shows that g = 1, the usual 
choice in the literature, is not a good scheme, but rather g = cxd is a good scheme to choose. 

Combining rriBR and /br in Eqs. (fT7|) and (jH^ . respectively, we finally obtain the B 
decay rate given in Eq. (pH) with 



To conclude, we summarize our results as follows. 

• The leading renormalon in the pole mass can be resummed accurately by employing 
a bilocal expansion. The BR mass ttibr, which is the real part of the integral is 
unambiguous and converges rapidly, although it has the same divergent perturbative 
expansion as the pole mass and is not a short distance mass. The renormalon ambi- 
guities in the pole mass cancel when a consistent integration contour is taken for all 
Borel integrals. Defined as the pole mass, the BR mass gives a precise meaning to the 
relation between the pole mass and a high energy mass like the MS mass. 

• A leading renormalon in the effective field theories is no longer a huddle but a blessing 
for the perturbation theory. By providing a manageable constraint on the functional 
form of the Borel transform it allows one to speed up the convergence of perturbation 
theory. 

• The renormalon residue of the pole mass can be determined accurately. Using the 
renormalon cancellation one can determine the residue of the other perturbative ex- 
pansion that cancels the renormalon ambiguity in the pole mass to the same accuracy. 
In the case of the B decay, we used this information on the residue to estimate the 
NNLO coefficients /a in Eq. ((221) • 

• The renormalization scheme dependence of /(cts), parametrized by Eq. ()25|) . is signif- 
icant and should be taken into account. 




(34) 
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• The B decay rate was calculated to NNLO (using an estimate of /2) in the BR scheme. 
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